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An electromagnetic (EM) Bloch wave propagating in a photonic crystal (PC) is characterized by 
the immittance (impedance and admittance) of the wave. The immittance is used to investigate 
transmission and reflection at a surface or an interface of the PC. In particular, the general properties 
of immittance are useful for clarifying the wave propagation characteristics. We give a general proof 
that the immittance of EM Bloch waves on a plane in infinite one- and two-dimensional (2D) PCs 
is real when the plane is a reflection plane of the PC and the Bloch wavevector is perpendicular to 
the plane. We also show that the pure-real feature of immittance on a reflection plane for an infinite 
three-dimensional PC is good approximation based on the numerical calculations. The analytical 
proof indicates that the method used for immittance matching is extremely simplified since only 
the real part of the immittance function is needed for analysis without numerical verification. As 
an application of the proof, we describe a method based on immittance matching for qualitatively 
evaluating the reflection at the surface of a semi-infinite 2D PC, at the interface between a semi- 
infinite slab waveguide (WG) and a semi-infinite 2D PC line-defect WG, and at the interface between 
a semi-infinite channel WG and a semi-infinite 2D PC slab line-defect WG. 



PACS numbers: 42.70.Qs, 42.79.Gn 

I. INTRODUCTION 

Photonic crystals (PCs) are artificial materials whose 
dielectric functions periodically vaiy±2£££iL The pri- 
mary characteristics of PCs are photonic band gaps and 
peculiar dispersion relations, both of which enable the 
propagation of electromagnetic (EM) waves in PCs to 
be controlled. Recent interest in PCs has focused on 
"open-system pCs,"^£iL£*2*i2iiI which are connected to 
other structures at their boundaries. Typically investi- 
gated are a finite PC contacting another medium and 
finite PC line-defect waveguides (WGs). In open-system 
PCs, the energy of propagating EM waves is transferred 
into or out of the PCs through the boundaries. This is in 
contrast to closed-system PCs, which are infinite in size 
and therefore contain all the energy at any time. 

An open-system PC can be divided into a set of sub- 
systems — the PC and the attached structures. There- 
fore, it should greatly facilitate designing open-system 
PCs if each subsystem is independently designed and 
then all the subsystems are assembled. Such a method 
has become well established for electrical circuit design. 
The behaviors of the electrical energy flows are pre- 
dicted based on electrical immittance (impedance and 
admittance)^ Two electrical circuits can be assembled 
without energy reflection at their interface if the electrical 
immittances of both circuits match. A complete electri- 
cal circuit system can thus be designed after designing 
each subsystem independently. 

EM wave media including open-system PCs can be re- 
garded as EM wave circuits, for which an analogous figure 
of merit, EM wave immittance, is similarly defined^ It 
is a ratio of the electric field to the magnetic field of a 
propagating or decaying wave in a single medium. In 
practice, it has been possible to design one-dimensional 



(ID) open-system PCs by using the EM wave immit- 
tance. Particularly notable is that the renormalized Fres- 
nel coefficient for a semi-infinite ID PC with an arbitrary 
unit cell has been defined using the EM wave immittance 
for Bloch wavesjM However, there is still difficulty in de- 
signing multi-dimensional open-system PCs. This is be- 
cause the number of EM wave immittances that should 
be considered at a point on an arbitrarily shaped bound- 
ary increases with the number of structural dimensions. 
Specifically, the immittance of an EM wave with a given 
polarization is expressed as one parameter for a ID PC, 
two parameters for a two-dimensional (2D) PC, and six 
parameters for a three-dimensional (3D) PC. Moreover, 
each parameter of an EM wave immittance is complex in 
general. Therefore, excessive computational cost would 
be needed to match every immittance parameter between 
two multi-dimensional subsystems. 

One way to simplify the matching problem is to re- 
duce the number of parameters by considering structural 
symmetry. Also, simplifying the parameters themselves 
is helpful, as has been suggested by Boscolo et al>i^ They 
pointed out that a propagating wave in a PC line-defect 
WG has a pure-real immittance on periodic symmetry 
planes normal to the direction of propagation. Similar 
findings have been obtained through analytical calcula- 
tion of the impedance of ID PCs^i^ and through nu- 
merical calculation for 2D PCs>i£ In this article, we take 
another approach to simplifying the matching problems 
— qualitative optimization of individual subsystems. 

This paper is organized as follows. In Sec. [H] we 
first prove that pure-real immittance is a common prop- 
erty of EM Bloch waves in ID and 2D PCs, hence it 
can be used for any ID or 2D open-system PC. We also 
prove that the pure-real property of the immittance is 
true for 3D periodic structures, but under a limiting con- 
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dition. Then, in Sec. IIHI we present three examples 
of designing multi-dimensional open-system PCs. This 
analysis can be made because the immittance defined by 
using eigenmodes on the boundary of two semi-infinite 
regions (a PC and its external structure) can be deter- 
mined from the eigenmodes of an infinite PC and its ex- 

irrng 



ternal structure. In Sec. IIII Al we present a method for 
qualitatively evaluating the reflection at the surface of a 
semi-infinite 2D PC; the approximate availability of the 
perfect antireflection coating (ARC) method developed 
for a semi-infinite ID PC— to a semi-infinite 2D PC is 
emphasized. These semi-infinite 2D PCs can be used 
for analyzing various applications such as add/drop mul- 
tiplexers, dispersion compensators, polarization filters, 
and image processors Ai&iii In Sec. IIII Bl we investigate a 
2D line-defect WG that efficiently transmits EM energy 
to a channel WG. The semi-infinite 2D PC line-defect 
WG is a promising platform for ultra-compact photonic 
circuitsASI I n Sec. IIII CI we investigate the interface 
between a semi-infinite channel WG and a semi-infinite 
2D PC slab line-defect WG as a typical 3D problem. 21 
All the applications of the proof are presented using the 
concept of immittance matching. Section llVI devotes for 
conclusion. 



II. GENERAL PROOF FOR IMMITTANCE OF 
ELECTROMAGNETIC BLOCH WAVES 

We begin our description of the proof with a general 
expression of EM Bloch waves. Assuming time harmonic 
EM waves, the general expression for an EM Bloch wave^ 

E(r,t) = u k „(r) exp(ik • r - iu kn t) , (1) 
H(r,t) = v k „(r) exp(ik • r - iu> kn t) , (2) 

where n and k stand for the band index and Bloch 
wavevector, and functions Uk n (r) and Vk„(r) have lattice 
periodicity. The immittance of a Bloch wave is closely 
related to the energy flow conveyed by the wave; it is 
therefore useful to define a complex Poynting vector: 

S(r) = iE(r,i) x H*(r,t) = iu kn (r) x v kn (r) , (3) 

where the real (imaginary) part of S is the density of the 
time-averaged (reactive) power flow— 

Here, we consider the time reversal state of Eqs. (QJ 
and J2J, which is written as 

t' = -t and (4) 

E( tr )(r,0 = E(r,-t) , (5) 

U^(r,t') = -H(r,-t) , (6) 

due to the even (odd) parity of electric (magnetic) field 
under time reversal operation^ Eqs. can be re- 

duced to£ 
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FIG. 1: Schematic illustration of typical infinite 2D PCs with 
reflection symmetry a y (at y = 0). (a) square lattice, (b) 
triangular lattice of circles with permittivity £2 and lattice 
constant a in medium with ei. 



Eqs. 0-© express the time reversal degeneracy of the 
spin 1 particle (photon). Eqs. Q-© will be used later. 

Next, we consider a Poynting vector described using 
EM Bloch waves under a reflection operation in infinite 
2D PCs. FigureQ]illustrates typical infinite 2D PCs with 
reflection plane a y . When the propagation direction of 
a Bloch wave is perpendicular to this plane, the electric 
Bloch wave E(r) is transformed by reflection operation 
into 



6 a E(r) = IC 2y u nky ] y (x,-y,z)e- lk yy, (10) 
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where Oy , /, and C 2y stand for the a y reflection oper- 
ator, the inversion operator, and the C 2 rotation about 
the y axis. Note that the obvious factor exp(— io>kni) is 
omitted in Eqs. I|10(l and Ijlip. Note also that this reflec- 
tion operation (a y ) is not the symmetry operation of the 
k-group at k y e y , but of the point group in an infinite 2D 
PC. Similarly, magnetic Bloch wave H(r) is transformed 
into 
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nk y ,y 



(x,-y,z)e^y. (12) 



In the derivation of Eqs. and iJUJl, we use Eqs. 

© and The non-zero components of EM waves in 
the Cartesian coordinate system are (H x , H y ,E z ) for TM 
and (E x ,E y ,H z ) for TE polarization— 

Accordingly, the Poynting vector described using the 
EM Bloch wave under reflection operation can be written 
as 



1 nk y ,x V nk y 



c /2 (for TM) 

s /2 (for TE) . 



(13) 
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On reflection plane y = 0, S y ay ^ (x, 0, z) equals S y (x, 0, z). 
Therefore, using Eqs. © and 113(1 . we obtain 



S y (x,0,z) = S^\x,0,z) 

^nk y ,x 
Vnky ,z 

lm(Z zx ) = (for TM) 
lm(Z xz ) = (for TE) , 



(14) 




where Z zx = u n k y , z /v n k y , x , and Z xz = 
(x,0,z). 



(for TM) 
(for TE) 

~^nk y ,x i Vnky 



(15) 
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Eq. ((16(1 means that the immittance is real in infinite 
ID and 2D PCs on a plane when the plane is a reflection 
plane and the Bloch wavevector is perpendicular to the 
plane. It also means that the reactive power flow is zero 
for the Bloch wave on the plane under consideration. 

For the case of an infinite 3D PC, Eq. ((14(1 can be 
written as 

Unky ,z^nk y .x ^nk y ,xVnk y ,z — 

(u nky , z v nky , x - u nky<x v nky , z y .(17) 

Therefore, an additional reflection plane is required for 
removing the imaginary part of the impedance. When 
there is reflection plane er 2 at z = 0, some components of 
the EM waves vanish on the plane. They are called TM- 
or TE-like based on the analogy to infinite 2D PCs^ In 
this case, the impedance is real on the intersection line 
of the y = and z = planes. 



III. APPLICATIONS OF PURE-REAL 
PROPERTY OF IMMITTANCE 

As a typical application of the above proof, we present 
a method reducing a reflection loss at the surface of a 
semi-infinite 2D PC, a semi-infinite 2D PC line-defect 
WG, and a semi-infinite 2D PC slab line-defect WG 
based on the concept of the immittance matching. The 
method for the immittance matching is extremely sim- 
plified since only the real part of immittance function for 
an infinite PC is needed for the analysis. The reflection 
coefficient of a semi-infinite PC should be determined via 
Bloch wave expansion in the semi-infinite PC and the ex- 
pansion should be described by all the eigen modes of an 
infinite PC (i.e., including decaying waves) ~ However, 
here we investigate only the impedance of propagation 
modes. This is another merit of our method, that is , 
it is not necessary to know all the eigenmodes in our 
method for immittance matching problem. The infor- 
mation on the impedance of propagation modes can be 
used for qualitatively investigating the reflection in vari- 
ous types of open system PCs. 
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FIG. 2: (a) Dispersion relation of FM wave in TM direction 
for infinite 2D PC with triangular lattice of air holes (t2 = 1). 
Permittivity ei is 11.9025 and hole radius is 0.4335a, where 
a is lattice constant. Index "A" ("B") indicates the coupled 
(uncoupled) modei— *— ■ *— ■ (b) Normalized impedance of Bloch 
waves with k y = 2PM/3 at y — 0. We use index B1-B5 (band 
number) to label the modes in order of increasing frequency. 
Only Bloch waves that coupled with an external plane wave 
are shown. 



A. Semi-infinite 2D PC 

The geometry of the first example is shown in Fig. ^ 
(b). The parameters of the PC is summarized in the 
caption of Fig. [21 

Figure [21 (a) illustrates the dispersion relation of the 
TM polarization wave for the infinite 2D PC in the TM 
directioni^i Suppose that one considers the reflection at 
the surface of a semi-infinite PC made from the infinite 
PC cleaved at the plane y = shown in Fig. [2 (b). The 
incident wave is a normal incident plane wave propagat- 
ing in vacuum. Then, one should calculate the immit- 
tance of the Bloch wave at y = by using the eigen mode 
of the infinite PC. The immittance of the Bloch wave is 
generally complex and has lattice periodicity of the in- 
finite 2D PC. The impedance at y — calculated using 
propagating Bloch waves with k y — 2rM/3 is shown in 
Fig. 121(b). The illustrated region is indicated by the thick 
solid line in Fig^ (b). The impedance is normalized by 
the impedance in vacuum. Note that only the impedance 
of Bloch waves able to couple with an incident plane wave 
("coupled modes") indicated by "A" in Fig. |2{a) is plot- 
ted. Since the spatial modulation of bands Bl and B2 is 
flatter than that of bands B4 and B5 and the impedance 
is near the value in vacuum (i.e. impedance normalized 
by vacuum value is near one), the reflection loss of bands 
Bl and B2 at the frequency is lower than that of bands 
B4 and B5 for the normal incidence of a plane wave. In 
addition, the flatness of the impedance enables us to re- 
duce the reflection loss by using an antireflection coating 
(ARC) method developed for semi-infinite ID PCs>i^ In 
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FIG. 3: Dispersion relation of infinite slab waveguide (WG) 
and infinite PC WG. PC has triangular lattice of air holes, 
with radius r = a/3 in medium with = 7.6f76. Even (odd) 
mode under reflection operation a x bisecting WGs is shown 
as solid (broken) line. Circles show points equivalent to K =- 
0.30, -0.35, and -0.40, which are used to calculate impedances 
shown in Figs. flja) and (b) . 

this case, since the plane y = is a reflection plane, 
thickness d and refractive index ni of the ARC can be 
estimated using 

D ee _l_ = l,3,5,--- and (18) 
Ao/4n 2 

n 2 ~ y/mBB(Y zx ) /e c at x,y = 0, (19) 

where D is the normalized thickness of the ARC, Ao is 
a wavelength in vacuum, n\ is the refractive index of a 
semi-infinite medium, and Y zx = 1 /Z zx for TM polariza- 
tion. Accordingly, the reflection loss is easily removed by 
using Eqs. Ijl8|l an d (|19|) . If a spatially flat impedance 
appears in the TE polarization case, Eq. (|19|) can be 
used by replacing Y zx with Y xz = 1/Z XZ to reduce the 
reflection loss of the TE polarization waves. 

B. Semi-infinite 2D Line-defect PC Waveguide 

The second example is related to the interconnec- 
tion between a semi-infinite slab waveguide (WG) and 
a semi-infinite 2D PC line-defect WG, or "PC WG." 
The connection between different WGs is important to 
realize the ultra-minute photonic circuits. Although 
some numerical investigation were made for similar 
structures) 20 i 21 i 31 i 32 i 33 i 34 i 35 we investigate it based on the 
concept of immittance matching. The semi-infinite PC 
WG is made by removing a row of air holes in the TK 
direction from a semi-infinite 2D triangular lattice of air 
holes with radius r — a/3 and e% = 1 in a medium with 
ei = 7.6176. The slab WG has width y/3a - 2r and 
permittivity e = 7.6176 (the same as the t\ of the PC). 
The dispersion relation of the propagation modes of the 
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FIG. 4: Dependence of spatial distribution of impedance on 
wavevector and on surface of semi-infinite PC WG. Differ- 
ence between Figs, fl] (a) and (b) is the position of the sur- 
face of semi-infinite PC line-defect WG connected with semi- 
infinite slab WG, shown as dotted line in diagram— above 
each graph. Wavevector dependence is shown as solid (nor- 
malized wavevector K = —0.30), dashed (K = —0.35), and 
dotted [K = —0.40) lines. The rectangular impedance distri- 
bution is impedances for slab WG with width \/3a — 2r. 



infinite slab WG and infinite PC WG are shown in Fig. 

The polarization is TE; that is, the magnetic wave is 
parallel to the slab and to the axis of air holes in the infi- 
nite 2D PCH Since the both of the slab WG and the PC 
WG have reflection plane bisecting these WGs, the prop- 
agation mode has even (broken line) or odd (solid line) 
parity under reflection operation, as shown in Fig. [3] 
The shaded region indicates the projection band struc- 
ture. When K= 0.6 ~ 0.9, the dispersion relation of 
the infinite slab WG is similar to that of the infinite PC 
WG. Such situation was called k-matching (wavevector 
matching) by Mekis et ali^^ 5 as a mechanism to reduce 
the reflection from the joint of a distributed Bragg reflec- 
tor WG and a 2D PC WG with square lattice of rods. On 
the other hand, the exactly derived renormalized Fresnel 
coefficients of a semi-infinite ID PC did not include Bloch 
wavevector explicitly^ The analytical investigation of a 
relation between immittance matching and k-matching 
will be published elsewhere^ 

Here, we investigate the direct connection between a 
semi-infinite slab WG and a semi-infinite PC WG based 
on the immittance matching. Figure 0] illustrates the de- 
pendence of the spatial distribution of the impedance on 
the wavevector and on the surface of the semi-infinite 
PC WG. The difference between Figs. H (a) and (b) is 
the position of the surface [shown as a dotted line in 
the diagram above each graph in Figs. Efa) and (b)] 
of the semi-infinite PC WG connected with the semi- 
infinite slab WG. For the semi-infinite slab WG, the 
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impedance is homogeneous inside and outside the slab 
region, which shows a rectangular spatial distribution. 
In addition, the impedance at a surface is independent 
of the cleaved position of the infinite slab WG. For the 
semi-infinite PC WG, the impedance depends strongly 
on the position of the surface of the semi-infinite PC 
WG due to the spatial periodicity of the impedance 
in the infinite PC WG. Moreover, the impedance on a 
surface has a complex spatial distribution on the sur- 
face. Divergence appears where the magnetic wave van- 
ishes. Divergence in a differential impedance appears 
at normalized position X = 2x/\/3a — ±2 ± 2/3a/3, 
±4 ± 2/3v / 3- • ■ for Fig. H (a), and at ±1 ± 2/3^3, 
±3 ± 2/3V3 • • • for Fig. H (b). The divergence is due 
to the discontinuity of the electric wave perpendicular 
to the air holes at the surface of the semi-infinite PC 
WG. From the Poynting vector calculation, we can con- 
firm that the EM energy of these propagation modes 
concentrates around the normalized position X = Oi^ 
Therefore, region X = -1 + 2/3\/3 ~ 1 - 2/3\/3 (-the 
width of the slab WG) is quite important for analyzing 
the immittance matching between the semi-infinite slab 
WG and semi-infinite PC WG modes. In this region, 
the impedance matching between the semi-infinite slab 
WG and the semi-infinite PC WG modes of Fig. El (b) 
is better than that of Fig. 01(a). Hence, the geometry of 
Fig. E] (b) is better for reducing reflection loss. In addi- 
tion, when the wavevector approaches the 1st Brillouin 
zone boundary, the difference in the impedance between 
the slab WG and the PC WG modes becomes notice- 
able. Therefore, in a direct interconnection between the 
semi-infinite slab WG and the semi-infinite PC WG, the 
reflection loss of modes with K = —0.4 is higher than 
that with K = -0.3. 



Application to 3D structures —Semi-infinite 2D 
PC Slab Line-defect WG- 



In this section, we investigate the connection between 
a semi-infinite 2D PC slab line-defect WG and a semi- 
infinite channel WG based on the concept of immittance 
matching. Schematic illustrations of two kinds of connec- 
tions are shown in Figs. El (a) and (b). The structures 
are sandwiched by two semi-infinite vacuum regions in 
the vertical direction (s-axis direction), hence they are 
3D in nature. The 2D PC slab consists of a triangular 
lattice of air holes. The line defect WG is introduced into 
the slab by filling the air holes in a row. The detailed pa- 
rameters of these structures are given in the caption of 
Fig. [SI The width and thickness of the channel WG are 
•\/3a — 2r and 0.6a, where r is the radius of the air holes 
and a is the lattice constant. 

Figure El illustrates the energy dispersion relation of 
the TE-like modes of an infinitely long 2D PC line-defect 
WG. Fundamental EM mode propagating in the channel 
WG can be coupled only with the "odd" modes of 2D PC 
slab WG in the ideal case [Mode symmetry is defined by 
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FIG. 5: Dispersion relation of TE-like modes of 2D PC slab 
line-defect WG. The slab consists of the triangular lattice of 
air holes (radius r = 0.29a, thickness h = 0.6a, refractive 
index n — 3.4, where a is the lattice constant). The line defect 
is introduced by filling up a row of air holes. These parameters 
are the same as in Ref. 133 . The 3D calculations are needed 
for this structurei— The symmetry of the non-leaky guided 



modes is indicated as cry-even or 



odd, where a y indicates the 



symmetry plane at y = in Fig. 03(a) and (b). Squares stand 
for points executing impedance calculations shown in Figs. [5] 
(a) and (b). 



the symmetry operation associated with the symmetry 
plane (a y at y=0) bisecting the 2D PC slab WG and the 
channel WG]. Therefore, we consider the odd mode only 
hereafter. 

Even if the structures of the 2D PC slab WG and the 
channel WG were fixed, there would still be ambiguity 
in determining the structure of the connection — the po- 
sition of the cutting surface in a unit cell of a 2D PC 
slab WG. This situation is the same as the 2D problem 
in Sec. IIII Bl If we select a reflection plane as the cut- 
ting surface, there are two possible choices, as illustrated 
in Figs. El (a) and (b) (at z = +0). To illustrate their 
difference more precisely, we show the top views (x = 
plane) and the side views (z = plane) of the refractive 
index distribution in the unit cell used for the supercell 
calculations^ in Figs. El (c)-(f); black indicates the re- 
gions with refractive index n = 3.4, and white indicates 
the region with vacuum. The size of the supercell used 
in our calculation can be obtained from these figures. 

When a 2D PC slab WG is terminated and connected 
with a semi-infinite channel WG at z = 0, the immit- 
tance distribution should be calculated at z = +0. The 
analytical results presented in Sec. HTlshow that the pure- 
real feature of immittance on a symmetry plane for a 3D 
structure is preserved along a particular line only. How- 
ever, our numerical results indicate that the imaginary 
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FIG. 6: Dependence of spatial distribution of impedance on wavevector and on surface of semi-infinite 2D PC slab line-defect 
WG. Difference between Figs. |5] (a) and (b) is the position of the surface of semi-infinite PC slab line-defect WG connected 
with semi-infinite channel WG (z — +0). (c)-(f): top and side views of refractive index distribution, (g)-(l): Dependence of 
impedance RefZj^] (in units of 1/eoc) on cutting surface and Bloch wave vector. The dotted lines indicate the width of the 
channel WG. 



part of immittance on a symmetry plane is much smaller 
than the real part. In this situation, the imaginary part 
of the immittance on a symmetry plane is negligible in 
qualitative immittance matching analysis. Moreover, ad- 
ditional components of the immittance appear in the 3D 
problem. This is because all components of an EM field 
cannot be treated separately. The two possible combina- 
tion of EM fields in the impedance calculation for Bloch 
waves propagating in the z direction are 



Z yx (x,y, +0) 
Z xy (x,y, +0) 



Hk z ,x(x,y,+0) 
E kz , x {x,y, +0) 



(20) 
(21) 



xy 



IS 



H k^ y (x,y,+0) 

Since we use TE-like mode only in this analysis, Z. 
much smaller than Z yx , and this was confirmed by our 
numerical calculations. Therefore, in Figs. (g)-(l) we 
illustrate the real part of impedance distribution Z yx at 



z = +0. The regions of these plots are the same as those 
in Figs. (e) and (f). The value of Re^J is colored 
linearly. Note that the region |Re[Zj,J| > 2 is colored 
red or blue. Note also that the energy dispersion relation 
of the odd mode shown in Fig. [S] has negative group 
velocity, so the impedance is calculated by using time 
reversal states of EM fields. 

The results presented in Sec. 1111 Bl clearly show that 
the structure in Fig. (a) is more suitable for low-loss 
connections than that in Fig. (b). This is supported 
by 3D calculations of the impedance distribution around 
the WG region < y < V3a — r (between dot- 

ted lines in FigEJ] . The impedance distribution shown in 
Fig. IHl(g) is flatter than that in (h) at this region; mean- 
ing that the structure in (a) has lower reflection than that 
in (b). 

Next, we investigate their dependence on the Bloch 
wavevector, which is also shown in Figs. El(g)-(1)- When 
the length of the Bloch wavevector increases, the regions 
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colored red and blue expand. This means that the re- 
gions with |Re[iL z ]| > 2 expand, so the reflectance of 
EM waves increases. 



IV. CONCLUSION 

In conclusion, we presented a qualitative method 
for evaluating the reflection of multi-dimensional semi- 
infinite PCs that is based on the concept of immittance 
matching. Using the analytical investigation of Poynt- 
ing vectors defined using Bloch waves, we presented a 
general proof that the imaginary part of the complex im- 
mittance on periodic planes is zero in infinite ID and 2D 
PCs when the plane is a reflection plane and the Bloch 
wavevector is perpendicular to the plane. The infinite 
3D PC case was also analyzed. To show the usefulness 
of this proof, we presented a qualitative method for eval- 
uating reflection loss (a) at the surface of semi-infinite 
2D PCs, (b) at the interface between a semi-infinite slab 
WG and a semi-infinite PC line-defect WG, and (c) at 
the interface between a channel WG and a semi-infinite 
2D PC slab line-defect WG. The first example showed ap- 
proximate applicability of the ARC method developed for 
semi-infinite ID PCs to semi-infinite multi-dimensional 
PCs. The second example showed that the use of a con- 
nection with the geometry shown in Fig. 0] (b) reduces 
the reflection. The third example showed that the pure- 
real feature of the immittance is a good approximation 
on a symmetry plane of 2D PC slab WG. 

Note that reflection of PCs is sometimes discussed 



in relation with the group velocity or the density of 
states in the infinite PC. While this can provide good 
physical intuition, we need to remember that the 
reflectance strongly depends on the structure of the 
surface of a PC. The immittance of EM Bloch wave is 
essential for the discussion of the reflection phenomena 
for any "open-system" PCs. The general concept of 
immittance matching was used for this theoretical 
prediction and is well suited for reflection and trans- 
mission analysis for a wide variety of "open-system" PCs. 
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